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We use an exact holon and spinon Landau-liquid functional which describes the 
holon - spinon interactions at all scattering orders, to study correlation functions 
of integrable multicomponent many-particle problems showing both linear and non- 
linear energy bands. We consider specific cases when the dominant non-linear band 
terms are quadratic and apply our results to the evaluation of the ID Hubbard model 
correlation functions beyond conformal-field theory. 
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Recent experimental investigations on the one-particle spectral properties of one- 
dimensional (ID) insulators |I[ have confirmed that the low-energy physics of these materials 
is dominated by holons and spinous and that the conceptual importance of these excitations 
is the same as that of other elementary particles and quasiparticles in solids such as phonons 
and magnons. The discovery of these elementary excitations followed the exact diagonaliza- 
tion of the ID Hubbard model H , whose Bethe-ansatz (BA) solution was obtained in terms 
of two coupled equations, associated with charge and spin modes. However, that solution 
does not provide direct information on matrix elements, and previous theoretical studies on 
the one-particle spectral properties of this model were mostly numerical |3|j4j]- Analytical 
studies on the problem either referred to the limit of on-site Coulomb interaction U — > oo, 
when the charge - spin separation occurs at all energies, or to the metallic phase away 
from half filling, where conformal-field theory (CFT) has provided important information 
on correlation functions [f|-§]. Unfortunately, that theory does not allow the evaluation of 
one-particle correlation functions for the half-filling Mott insulator and the study of these 
functions at finite values of U remains an open problem of great physical interest. 

In this Letter we combine a Landau-liquid functional which describes exactly the holon 
and spinon interactions at all scattering orders with new detected symmetries, to derive 
asymptotic expressions for the ID Hubbard model one-particle correlation functions both 
at half filling and maximum magnetization. Our critical theory refers to all integrable 
multicomponent quantum liquids whose low-energy physics is dominated by elementary 
excitations with both nonlinear and linear energy band terms. We study the case of quadratic 
bands in detail and also consider the general case of nonlinear band edges. Our results reveal 
surprising invariances of great importance for the theoretical understanding and description 
of the unusual spectral properties detected in the new low-dimensional materials. 

We consider the ID Hubbard model in a magnetic field H and chemical potential fi which 
describes N = iVj + N± interacting electrons and can be written as H = T + U[D - ±JV] + 
T, a =c,s Va^^z where f = -tJ2j,a[ c ]a c j+i° + h.c] is the "kinetic energy", D = J2j n j,^ n j,i 
measures the number of doubly occupied sites, ^ = /i, /i^ = fioH, /i is the Bohr magneton, 
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S c z = ~i[N a ~ N), S s z = -±[A> T - A\], N = £ CT A\, N a = Ejn j>(T , h ji(J = c] a c ja , c) a and c ja 
are electron operators of spin projection a at site j = 1, N a , and t is the transfer integral. 
We choose a density n = n-f + nj_ = N/N a and a spin density m = n-j- — rz| in the intervals 
< n < 1 and < m < n, respectively, where n a = N a /N a and consider kp a = irn a and 
hp = irn/2. We use units such that a = h = 1, where a is the lattice constant. 

It is convenient to use the operational representation for holon and spinon excitations 
introduced in Ref. [7J] which is exact and refers directly to the BA solution. In this repre- 
sentation the holons and spinons are the holes of the c and s anticommuting pseudoparticles 
||, respectively. It should be emphasized that this connection is only true in the limit of 
zero magnetic field and that the concept of c and s pseudoparticle is more general than that 
of holon and spinon and refers also to finite values of H. In this basis the Slater-determinant 
levels of a pseudoparticles, with a — c,s, are exact energy eigenstates of the many-electron 
problem and the total momentum is additive in the pseudoparticle momenta q. On the other 
hand, the Hamiltonian contains pseudoparticle interaction terms, yet integrability implies 
that all the interactions are of zero-momentum forward-scattering type. The exact many- 
electron ground state is in this representation described by the a-pseudoparticle distribution, 
KiQ) = - <?) and N®(q) = G(q - q^) with < q < q^ and q^ < q < 0, re- 

spectively, where = ±g^ a + 0(1/N a ) and qp a = irN a /N a . Here N c = N and N s = N± 
and hence qp c = 2kp and qp s = k F ±. The pseudo-Brillouin zones are g^ 1 ^ = ±q a + 0(1/N a ) 
and q a = irN*/N a where N* = N a and N* = iVj and hence q c = tt and q s = kp^. Note that 
at half filling, AT? = N* — N c = (zero magnetization, = N* — N s = 0) and the ground 
state has no holons (spinons). Importantly, N a and A^* are always conserving numbers. The 
eigenstate energy expressions only involves the momentum distribution N a (q) and for states 
differing from the ground state by a small density of excited a pseudoparticles the excitation- 
energy spectrum can be expanded in the momentum deviations, AN a (q) = N a (q) — N®(q), 
what leads to a Landau functional 0, AE L = YhIi AE iy which up to third order in i reads 

AE L = J2^N a (q)e a (q) + j £ £ AN a (q)AN a ,(q>)± f a>a ,(q, q') + 

a ,Q a,a' q,q' 
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+ 71 E E AN a (q)AN al (q')AN a //(q")-g aoi , a „(q,q',q") + fo.o. . (1) 

a, or, a" q,q ,q 

Importantly, the pseudoparticles do not map onto electrons as U is adiabatically turned off. 

All the coefficients of the functional ([I]) can be exactly derived from the BA equations 
||, yet for orders i > 2 the calculations become extremely lengthy. So far only the i = 1 
band e a (q) and i = 2 /q, jQ ,/ function were derived ||. ^From the band expressions the ve- 
locity v a (q) = de a (q)/dq and the function a a (q) = dv a (q)/dq follow. The f a ^ a i function 
expression only involves the velocity and the two-pseudoparticle phase shift Q a ,a' {q-, Q 1 ) @ • 
This phase shift plays a central role in our critical theory - it is an interaction-dependent 
parameter associated with the zero-momentum forward-scattering collision of the a and a 1 
pseudoparticles of momentum q and q', respectively. For instance, the phase-shift combina- 
tions £ l aa , = 5 a<a i +J2j=±i(jY®aa'(<lFa, jlFa 1 ) with % = 0, 1 fully determined the a anomalous 
dimensions of the generalized critical theory, as we discuss below. 

When both the densities of a pseudoparticles, n a , and of a pseudoholes, n^, are finite 
the dominant e a (q) band terms are linear in the momentum for low-energy excitations. We 
emphasize that in this case only the terms of order i — 1 and i = 2 of the functional (0) are 
relevant for the low-energy physics. Importantly, replacing in expression (H) the bands and / 
functions by the corresponding linearized bands and values of the / function at the pseudo- 
Fermi momenta, using a general deviation expression describing all low-energy processes, 
and performing the q summations leads to the CFT spectrum derived in Ref. by a finite- 
size analysis. This reveals that the pseudoparticle basis is the natural representation for 
the critical-point operator algebra. Moreover, for these densities CFT provides asymptotic 
correlation functions expressions ||. 

We emphasize that the functional (P is universal for a wide class of integrable multi- 
component models, with arbitrary number of a pseudoparticle branches. It occurs often in 
these systems that as the density n a or becomes small the dominant low-energy a band 
term is non linear and of the general form 

e a {q)n±\\q\-Q a \ j /j\M a (2) 
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where Q a = when n a — > and Q a = q a when — > 0, A^a is a positive generalized 
mass, and j is a positive integer, j = 1, 2, 3, 4, oo. Below we use the notation {a,j} 
for a a band of type j. In the present model there occur band changes {c, 1} — > {c, 2} as 
n c — > and nj? — > 0, {s, 1} — > {s, 2} as n s — > 0, and {s, 1} — > {s, oo} as Z7 — > oo (e s (g) = 
implies {s, oo}). A central point is that for each {a, j} branch the low-energy critical-theory 
spectrum has contributions from terms up to order % — j + 1 of the functional ([!]) . 

In this Letter we consider explicitly the low-energy critical theory associated with the 
case when some of the a branches have {a, 2} bands and for all the remaining branches 
(whose minimal number is one) the band is of {a, 1} type. This requires the use of terms up 
to the order % = 3 in the functional (|]). Lengthy calculations involving substitution of the 
third-order deviation expansions in the BA equations of Ref. || lead to g a , a ',a"(q,q',q") — 
9a,a',a"(q, <?', q") + 9a,a',a"(q', q\ q) + g(q", q, q') where 

g a ,a', a "{q, q', q") = 4(7r) 2 {a a (g)$ QjCt /(g, q')$ a ,a"{q, q") 

+ Y J2j-^-^a"',a{jqFa"', <?)$«'" ,a'{jqFa"> , q')$a">,a"(jqFa">, <?") 
j=±l a'" ^ 

a, a', a" (?> q', q")] 

+ Y V a"'^a"'AjqFa"',q)A a "' ta ', a "{jqFa"',q',q")} , (3) 

j=±l a'" 

A a ,a',a"{q,q',q") = ®a,a>(q,q')-j-®a, a "(q,q") 

dq 

, 9 / (o)/ n\ d n <s> a >Aq',q)<&a>,a»{q',q") ®a"Aq",q)$ a », a >{q",q') ] 

Pa iq)) dq i[ 2np* (r£V)) + 2 W (r^V)) 1 

j=±l a'" lTip a "> {r a '"j 

v a = v a (qF a ), a a = a a (qF a ), the function 2np a (r) is such that l/2np a (r^(q)j = 
drW(q)/dq, r^(q) = K<®(q), r®(q) = S<®(q), K^(q) and S^(q) are the rapidity functions 
obtained by solving the integral equations (Al) and (A2) of the second paper of Ref. @] for 
the ground-state distribution N°(q), and r a = r^\qp a ). 

The Landau character of the spectrum (|I[) implies that the only small parameter which 
defines the range of validity of the corresponding low-energy critical theory is the density 



of excited pseudoparticles. The associate low-energy processes can produce changes in the 
numbers N aL of right (a = +1) and left (a = —1) pseudoparticles (with i = sgn(q)l = ±1), 
which are conserved independently. Equivalently, they can produce changes in the charge 
and current numbers N a = E t N aL and J a = (1/2) E t '^m, respectively The above processes 
can be classified into three types, (i) ground-state - ground-state transitions associated with 
variations AN a = N a — N® which only change q Fa , (ii) finite momentum K = ^ a T> a 2qF a 
processes associated with variations V a = J a — J°, and (iii) a number Ngfc = 0, 1,2, ... of 
elementary particle - hole processes around the Fermi point qp a . In the absence of transitions 
(i) we define qFa relatively to the initial ground state and otherwise relatively to the final 
ground state. Independently of the form of the band e a (q), at low energy the above processes 
lead to a momentum of the form K = k + E a AP a where AP a = j^- [AN a V a + E t lN*£] and 
in the present case, k = E a T^a^QFa- Let us denote by Em Ea> an d E a the summations over 
linear branches, quadratic (or other j > 2 non-linear) branches, and all types of branches, 
respectively. The {a, 2} band reads e a (q) = ±||g| — Q a | J /2m*, with M. a = m* and m* = 
l/a a (0) and m* = l/a a (q a ) for pseudoparticles and pseudoholes, respectively. These masses 
play the same role as the effective mass of Ref. ||. Introducing in Eqs. (|l])-(f§) the {a,j} 
band expressions (with j = 1 or j = 2) and the low-energy deviations associated with all 
low-energy processes and performing the q summations leads to the critical spectrum 

As = ^LX>« A a + E«-ra + E{(^) 3 ^«) 2 + (^) 2 ^ + T^> ' ( 5 ) 

where Af = E i= o,i(A^) 2 + Nf , N» = E a '£° Q ^, K = Ea'CL^a', v a = v a C with 
v a = qFa/m a and v a = q Fa /m* a (or v a = [q a - q F a}/m a and v a = [q a - qFa}/m* a ), £° = 
^ = ^ = and £ a is independent of v a and contains no pseudoparticle interaction 

parameters. Note that qp a = 7iAN a /N a (or q a — q Fa = —7iAN a /N a ) in the limit when 
n a = (or n^=0) for the initial ground state. 

When all a bands are linear, the critical-theory expressions involve the anomalous di- 
mensions, 2A^ = £ i=0) i(.A^) 2 + iAN a V a + 2Ng L (equal to (E<=o,i ^a) 2 + 2iV 5 when there 
is dual symmetry) and the asymptotics of correlation functions is given by CFT || and reads 
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X&(x,t) oc Yi a L [e~ lk ° x /(x — iv a t) 2A ' a }. On the other hand, the critical theory associated with 
the low-energy spectrum (|5]) is characterized by two remarkable properties: (I) in spite of 
its non-linear velocity terms, the same above anomalous dimensions are obtained from the 
following equation, 2A L a = £-\J2a dAE/dv a + J2& dAE/dv a + iAP a ]. Since their expressions 
are the same as in the linear regime, their values are given by taking the corresponding 
limits of the densities in the general expressions - in the present model the values of the 
dimensions 2A L C (or 2A^) are in the quadratic case given by the limit n — > 1 (or m — > n) of 
the corresponding dimensions of the linear regime; and (II) as n a — > or n 1 ^ — > there is an 
adiabatic symmetry JTOj] which implies that , — > 5 a>a i in these limits. (However, if 
the a' band is linear, and remain in general interaction dependent.) This symmetry 
implies that £° — > 1, i.e. the pseudoparticle interactions which renormalize the quadratic 
mass m* and associate velocity v a vanish (and m a = m*). For instance, in the present 
model £° Q , ^ c -> 5 c ,q as n -> 1 {i.e., n h c -> 0) and £® a , ^ s -> S SjCt as m -> n (i.e., n s -> 0). 

These symmetries are observed in all multicomponent integrable models and imply that 
the terms of the critical spectrum (|5|) which contain no linear velocity terms are of non- 
interacting pseudoparticle character and only the linear terms in the velocities are inter- 
acting. Following this remarkable property and associate above symmetries, we find the 
following general expressions for the asymptotics of correlation functions 

X*(x,t)<x [I e-*»7[(x-M 2A «'^) E ' 2A "], x/pt/ml » 1, ae& 

a, a', i,' 

oc [J e- ik0X /[(-L'v a d) 2AL «> ( V^/m*)^ 2 ^] , x/t = 0. (6) 

These x and t dependences can be understood in the following way. In the limit of low energy 
each a excitation branch corresponds to an independent momentum-energy tensor compo- 
nent and to one independent Minkowski space with light velocity v a . For bands with both 
finite n a and densities the velocity v a is also finite and the metric is Lorentzian. How- 
ever, for vanishing small values of n a or the metric becomes Galilean. In the Lorentzian 
case the asymptotic of correlation functions involves the variables (x ± v a t) associated with 
Lorentzian transformations. On the other hand, in the case of Galilean symmetry space 
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x and time t are transformed independently. This is consistent with the two asymptotic- 
expression regimes involving either x or t. Moreover, the only combination of the time t and 
mass m* with dimensions of x is const x w't/m* . Note that the changes in the asymptotics 
only concern the metric whereas the a anomalous dimensions, whose values depend on U/t, 
n, and m, remain the same. Both the CFT asymptotic correlation function expressions and 
expressions (f)|) are limiting cases valid for finite and vanishing, respectively, values of n a or 
n\. As n a or n h a is gently increased, we come into a small-density {a, 2} — > {a, 1} band 
transition regime which is not described by these asymptotic expressions. 

Our results can be generalized for low-energy j = 1 and j > 1 bands of general form (pj). 
However, our theory does not describe the case when all bands are of j > 1 type. Although 
for j > 2 we have not calculated explicitly the coefficients of the functional ([!]) up to order 
i = j + 1, symmetry considerations lead to expression (H) for x/(t/ Aia) 1 ^ >> 1, whereas 
the x/t = expression becomes x&(%)t) oc Ila.av e~ tk ° x /[(— i'v^t) 2 ^ 1 (t/M a ) ~ j ]• Here 
j is meant to be a function of a, i.e. different a bands may have different j > 1 values. 
As we refer below, these expressions lead for the Hubbard model U — ► oo {c, 1} and {s, oo} 
bands to the same correlation function expressions as Ref. ||. 

Fourier transforms of the above asymptotic expansions provide correlation function ex- 
pressions for values of momentum k close to ko and low values of energy u measured from 
the initial ground-state energy. For instance, let A a (k,u), N a (k) = XL Ar(fc, and 
D a {uS) = J2k A$(k, uj) be the one a electron (or hole) spectral function, momentum dis- 
tribution, and density of states of the Hubbard chain. There are no previous analytical 
results for these functions at finite values of U and for (I) the half-filling Mott-Hubbard 
insulator (n = 1) and (II) the fully polarized ferromagnetic (rij = and = n = m) 
initial ground states. In case I the bands are of {c, 2} and {s, 1} type. For m = and 
both n = 1 and small finite densities of holes 5 and low negative (positive) values of uj 
and/or values of k close to kp our generalized theory leads to to A a {kp, oj) oc |co>|~ 7 / 8 and 
N a (k) oc | A: — A;^! 1 / 8 for the particle (hole) spectral function and momentum distribution. 
(At n=lwe define the ground-state energy at the bottom and the top of the Mott-Hubbard 
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gap for particles and holes, respectively.) On the other hand, the particle (hole) density 
of states is given by D a (u) oc | | 3// 16 and D a {uj) oc \uj\ 1 I s for n = 1 and small finite 6 
values, respectively. However, the latter expression (also predicted by CFT) is restricted to 
frequencies |u;| < E c = j(j) 2 32 5 m , . In the limit of n — > 1 this domain shrinks to a single 
point and the spectral function diverges as D(u) oc | | 3/ie _ j n cage jj bands are of 
{c, 1} and {s,2} type. We consider creation of one j electron and find for both m = n 
and for small finite densities n±, A±(kFi,w) oc | | — 1+ i C 1 — ^o] 2 anc [ N±(k) oc \k — fcpjj 2^ 1 ~' r i°^ 2 ^ 
where rjo = (2/tt) tan _1 ([4t sin(7rn)]/[7) and n < 1, whereas the density of states is given 
by Di{uj) oc |co>| - 5+5[ 1 - ? w] 2 anc [ D^(u) oc [cc?| 2 1 1-7 ? ] 2 for m = n and small finite values of nj,, 
respectively. The latter expression is restricted to energies \uj\ < E s = ^(x) 2 ^- Again, 
in the limit of —* this domain shrinks to a point and the spectral function behaves as 
Di(u) oc |c^| — i^i^ 1 ^o] 2 . These D±(u)) expressions are not valid for U — > 00 because then the 
bands are of {c, 1} and {s, 00} type and instead D^(u) oc While for n < 1 all our 

expressions lead to the correct results as U — > 0, the pseudoparticle theory also describes 
correctly the n = 1 Mott-Hubbard transition at U = 0, the c band being of {c, 2} and 
{c, 1} type for U > and at U = 0, respectively. Finally, concerning the comparison of 
our results with previously obtained m = and U — > 00 expressions, while our theory 
does not apply to the n = 1 case [fj of bands {c, 2} and {s, 00}, it provides expressions for 
the {c, 1} and {s, 00} case which corresponds to finite values of n and 5. Importantly, our 
general expressions lead to the same results as Ref. |4[], with the density of states given by 
D a {uj) oc |c<j|~ 3 / 8 and D a (uj) oc Ito] 1 ^ 8 for U — > 00 and small finite values of 4t/U, respectively. 
Our results reveal that the same anomalous dimensions control the low-energy critical theory 
independently of the linear or nonlinear character of its elementary-excitation energy bands 
and are expected to shed new light on the unusual properties of quasi-lD materials. 
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